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The possibility of ignition of some particles of a fragmented Under-Water (U/W) 
Shaped-Charge (S/C) jet is examined. It was determined that only minute particles (less 
than 6 micron diameter ) are likely to undergo complete melting. This estimate was based 
on the high stagnation temperature upstream of a particle moving at a speed of 7000 m/sec 
in water vapor within the cavity. 

Some two-color pyrometry tests of exploding wire that had been conducted at the 
Naval Postgraduate School were re-examined. It is suggested that color photographs of the 
corona region exhibit on going gas phase combustion. Due to the high temperature 
(inferred from the color), it is contended droplets do not occur in that region (either metal 
or oxide). This mode of combustion is distinctily different from metal droplet combustion 
such as may take place in U/W S/C or in rocket combustion. 

Existing models for metal droplet combustion were found to lack two features, namely 
full time-dependence and oxide cloud radiation for a large volatile droplet. We suggest to 
incorporate these features in the proposed model. 

For simplicity, only the one dimensional ( spherically symmetry ) case is considered. 
Metal liquid to vapor interface is not treated explicitly but through including the vapor qual- 
ity and vaporization latent heat into the total metal enthalpy. The condensed oxide (accord- 
ing to local equilibrium composition), is considered as a second phase. Thus the model 
ought to be adequate for either volatile or non-volatile process of metal combustion. 
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1. INTRODUCTION 



Metal particle ignition and burning phenomena have been used for several purposes. 
In solid rockets metal ingredients enhance the specific impulse and stabilize acoustic 
waves. In under-water vehicles metal particle burning is used as fuel for its motor. 
Recently explosive enhancement of Under-Water (U/W) Shaped Charge (S/C) by means of 
aluminum particle burning in water has been examined, (Ref. 1). 

In the case of stand-off explosion of an U/W S/C, a diverging shock wave propagates 
into the surrounding water. A hypersonic metal particle jet is moving in a radial direction, 
from the exolosive center to the target. As a result of severe heat transfer conditions at 
the stagnation point ofthe S/C metal jet, some particles may ignite. Ignition will depend on 
initial particle size, environmental condition and time exposed to heating fluxes. Combus- 
tion will enhance the pressure driving the shock wave and tear the target vigorously. Vari- 
ous combustion parameters affect the impulse on the target in an U/W S/C. Some of them 
are : heat of combustion of metal with steam, initial size of burning particles, flow-field etc.. 
A close investigation of different type of metals burning in steam was carried out by 
Exploding Metal Wire (EMW) tests, (Refs. 2-4 ). The purpose of these tests was to find the 
best available metal, which due to its burning, will improve the U/W S/C effectiveness. 
Although, as shown later, there is no similarity in combustion behavior between the U/W 
S/C and EMW, one can assume that the most effective and energetic metal found in the 
EMW tests, will also be the best to use in U/W S/C. 

Ref. 1 describes some experiments with aluminum shaped charge and a calculation of 
impulse improving by factor of 1.2 due to complete droplet combustion. Strong transient 
effects predominate the S/C particles heating process; therefore, ignition and combustion 
must be carefully examined. In some cases large particles would not ignite, (Ref. 5), and in 
others existing burning particles may be extinguished. The problem of predicting the 
amount of metal particle, varied by size, that contributes by combustion to the explosive 
enhancement, needs to be solved by constructing an appropriate model. Before postulating 
a new model, we need to understand the various parameters and stages in the general 
combustion process. We will discuss the metal ignition problem of U/W S/C in water vapor 
environment. Then we will examine the the governing parameters in the fast combustion 
process and construct an appropriate model by starting from the differential equations that 
describe the burning process. 
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2. GENERAL STAGES IN COMBUSTION 



There are three distinct time intervals in the combustion process: ignition, 
burning and cooling. The ignition stage includes heating of the metal surface to a certain 
level of temperature where slow chemical reaction will be replaced by near infinite reac- 
tion rates. As a result of rising wall temperature, surface reaction rate increases and thus 
the oxide layer grows. Heat is generated by the reaction, and the progress of the burning 
process depends on positive supply of heat to the metal. If heat loss rate is larger than the 
rate of heat released by combustion, the burning will be extinguished. The equation that 
describes the heat transfer balance is given in Ref. 5 as the following: 



(1./3.)p„c rdHd\ = - (k/2r)Nu(T-T ) 



( 2 - 1 ) 



Here 5 is the thickness of the oxide layer and can be presented by an Arrhenius type 
equation. Eq. 2-2 from Ref. 5 describes the kinetic law of oxidation. 



^5/5t = (B Exp(-E/(RT)) 



( 2 - 2 ) 



Spontaneous high heating rate, like an explosion or fast electrical heating, may cause igni- 
tion with a very small growth of oxide layer. Fig. 1 illustrates the difference in temperature 
rise between the fast and slow heating process. In the spontaneous case, fast heating rate 
raises temperature high above the ignition point and then, when the combustion is starting, 
the temperature drops below the level of adiabatic flame temperature. In contrast, com- 
bustion due to low heating will occur when the metal surface temperature reaches the 
ignition temperature. 



2 




boiling 

^^Xfteiting 

Mgnition 
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Fig. 1 Temperature rise in fast and low heating rates 
Some materials such as aluminum and magnesium have a Metal Boiling Temperature 
(MBT) lower than its Oxide Boiling Temperature (OBT). In such volatile material, combus- 
tion begins just after the metal surface is exposed to the surrounding atmosphere. Expo- 
sure of metal surface occur due to surface-tension-induced movement of liquid oxide layer 
gliding over the liquid metal or by blow-off of the oxide layer as happens in the exploding 
wire case. Vapor phase reaction will take place in the flame zone as long as heat losses 
are lower than heat generated by the combustion process. Classman stated in Ref. 6, that 
the maximum temperature in the combustion zone is bounded by the oxide evaporation 
temperature. 
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Non-volatile metals such zirconium or tantalum has MBT greater than OBT, and since 
the combustion temperature is lower than OBT, the metal would not be able to evaporate; 
therefore, combustion will occur on the metal surface, and the oxide product, near the 
reaction zone, will appear on the gas phase only. Table 1, from Ref. 7, presents some data 
on MBT and OBT for several useful metals. 



Metal 


Melting 


Boiling 


Oxide 


Melting 

Pt.,°C 


Boiling 
• Pt,,°C 


Adiabatic 
Flame Temp . 














(Equilibrium) 


Na (Sodium) 


. 97,8 


883 


Na20 


— 


1770 - 


i827 


Mg (Magnesium) 


650 


1103 


MgO 


2800 


3100 


3077 


Ca (Calcium) 


850 ' 


1690 


CaO 


2590 


3500 


3527 


A 1 (Aluminum) 


660 


2500 


Alj O3 


2050 


3500 


3527 


Ti (Titanium) 


1730 


3260 


TI 2 O3 


2100 


3000 


3027 


Zr (Zirconium) 


1845 


4900 


ZrO^ 

2 


2690 


4300 


4527 



TABLE 1 Physical properties of some metals and their oxides (Ref. 7) 



The composition of the environment has a strong influence on the combustion process. 
Static pressure as well as oxidizer concentration and type of ingredients determine the rate 
of combustion. Burning of an aluminum particle in free air is governed by vapor phase 
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reactions, whereas in rocket motors various species such as CO2 add some surface reac- 
tions to the vapor phase burning droplet (Ref. 8). Volatile metal burning is controlled mainly 
by vapor phase reactions and produce large quantities of small condensed oxide particles, 
often called smoke oxide cloud. The small condensed particles (less than 2 micron. Ref. 9) 
are produced in the flame zone and convected away by the main flow. If the flame zone is 
close enough to its moving burning droplet, smoke oxide particle may impinge on the 
metal surface, deposit on it and change its mass and momentum. For large particles the 
effect of smoke oxide radiation may not be neglected. King, in Ref. 8 neglected radiation 
effects and therefore limited his solution to a maximum of 30 micron aluminum particle 
size, and found good results compared with experimental tests in wet oxidizer. Ref. 16. 
Wood, Ref. 10 observed strong radiation effect on burning magnesium particles greater 
than 80 micron. Heat flux radiation, emitted from the condensed cloud, cools the reaction 
zone and heats the evaporating metal surface. A schematic description of the various 
zones in a burning moving aluminum particle, is given in Fig. 2, taken from Ref. 9. 




Fig. 2 General nature of the combustion zone of a volatile 
metal droplet burning in a mild convective flow ( Ref. 9) 
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3. ESTIMATING THE MAXIMUM PARTICLE SIZE FOR IGNITION IN AN U/W S/C 



The first step for estimating impulse enhancement of an U/W S/C by means of alumi- 
num combustion is to assure that the moving particles in the S/C jet do ignite under the 
special conditions. 

3.1 Conditions and Method of Estimation 

In order to make a quick estimation whether an aluminum particle can reach the burn- 
ing phase condition, while moving at high speed inside a shaped charge jet of particles 
cloud at an underwater explosion conditions, we evaluate the time needed by a single par- 
ticle to increase its temperature to the melting point and become completely molten. If 
that time, assuming severe conditions, is larger than the transit time from the stand-off 
explosion to the target then it is impossible for burning conditions. 

Time dependent conditions are not taken into account, and a simple calculation is 
being done by assuming an equilibrium state between water and steam, which means that 
the cavity is full with saturated steam. These conditions are the most severe boundary 
condition the shaped charge particles can meet. 

Other assumptions made here, are: 

-The surrounding pressure is equal to the hydrostatic water pressure. 

-All particles meet the same conditions and donot exchange heat betwe en themselves. 

-The particle spherical shape does not change. 

-Particle velocity remains constant and equal to 7000 m/sec. 

-Internal temperature distribution is ignored. 

3.2 Equations 

For the solid state particle the equation of lump heating is: 



X (T -T)-(j(€ T** - a T ^) = -p, 

c' r ' ' w g g I 



Ms'-p, 



Ms 



(3-1) 



Here X ,T ,T and a will be calculated according to the steam conditions and the particle 

c’ g’ r g ® 

velocity. The second equation shows the change in particle radius during melting time 
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Here >,^,7^,7^ and will be calculated according to the steam conditions and the particle 
velocity. 7he second equation shows the change in particle radius during melting time 



■ dydX (3-2) 

3.3 Calculation of the Coefficients 

a. Saturated Steam Condition at 3 atm. 

7 = 407 K ; p = 1.68 kgm/m 

R^= 848/18= 47.11 kgf*m/kgm* K 




7he speed of sound, a: 

a = = 771.32*47.11*9.81*407 = 498.38 m/sec. 

Viscosity at 7^: = 15.*10'” kgm/m*sec 

7hermal conductivity at 7^. k = 6.5*10'® kcal/m*sec* K 

Radiation absorbance, assuming equivalent length of 2 meter and a pressure of 3 atm: 

«g(3atm*2J= 0-7 

b. Recovery 7emperature and Stagnation Conditions. 

M = u/a = (7000. / 498.38 ) = 14. 

= T (1-t- (7-1) rlVlV 2 ) 

where r is the recovery factor. Assume r= 0.89, for turbulent flow. 

7j. =407. (1 + 0.16*0.89*14**)= 11766 K 
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This is a very high temperature, and we can, therefore, expect dissociation and an 
improvement in the thermal conductivity and viscosity. The density, viscosity and Prandtl 
number at the stagnation condition are: 

p° =( p/ (R^T^) = 3. 10^^/(47.11*1176.) = 0.054 kgm/m^ 

=40.*10’® kgm/m*sec 
Pr“ = 0.88 

c. Aluminum Properties. 

= 2700 kgm/m^ 

p^^ = 2380 kgm/m^ 

Pms ■ PMf " kgm/m"* 

= 932 K 

hj^^^ = 94.45 kcal/Kgm 

c = 0.23 kcal/kgm* k 
’’ms 

kg = 0.05 kcal/m*sec* k 

d. Stagnation Heat Transfer Coefficient. 

The heat transfer coefficient is defined as: 

X Pe Ue C St 

P« 

St is the Stanton number for stagnation point flow around a sphere 

(\5 .0^ 0 (V6 0 OS 0 05 

St = 0.763 (I5D/U ) (Re^)) (Pr )' (p /p) (p /p) 

For Mach number greater than 10, D|5/u = 1. 

In order estimate the shortest time for melting, we will assume a small particle diameter 
of the size of 20 micron, the appropriate Re^ number will then be: 

ReQ =1.68*7000*20’ 10'® /1 5* 10’^ =15680. 

The ratio of density at stagnation and free stream and the ratio of viscosities at these two 
conditions, are: 

p°/p =0.054/1.68= .032 ; p% =40/15= 2.67 
The heat transfer coefficient is; 

X^= 1.68*7000* .453*0.763*(15680)’(*'*)*0.88’‘^*(0.32*2.67)°^= 10.23 kcal/m'sec K 
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3.4 Comparison Between Forced Convection and Radiation Heating 

The convection term is defined as: 

qc = MTr-T) 

Assume particle initial temperature of 500 K 

= 10.23'( 11766- 500.)= 1.153*1o’ kcal/m^*sec 

Maximum radiation heat losses that can occur at the metal surface are at 

the highest temperature, T = T^ 

=<T(^^T^‘‘-agTg^) = 1.36*10"” (0.2*932“ .- 0.7*407“) = 1.46 kcal/m" *sec 
As q^ is much greater than q^. we may neglect radiation heat losses. 

3.5 Diameter Limit for Ignoring Temperature Distribution 

Temperature gradient inside matter can be neglected when k^/5i*r is greater than 6. In 
that case : 0.05/10.23*r is greater than 6. The condition is valid for every particle diameter 
less than 1630. micron. 

3.6 Calculation of Heating Time 

The calculation of the heating time needs the solution of the two differential equations 
(3-1) and (3-2) dropping the radiation term. Without the radiation form the solution is easy at 
two different periods of time. 

a. Conduction in the Solid State Time Interval. 

Equation (3-1) can be written for a non-dimensional variable O, as follows: 

(dO/dt) - t"'o = 0 (3-3) 

Where: O = (Tj^^- T.)/( T^-T.) ,T. is the initial temperature, and t is the time constant and it 
is equal to : p„ c v , 3*X 

^ po/ c 

The solution for equation (3-3) is: 

0=1- exp(-t/t) 

If tpp is defined as the heating time for a particle to reach the nondimensional melting 
temperature O^, we can find, from the solution of equation (3-3), the ratio of ® 

function of the volumeteric heat capacity Pj^^c and the heat transfer coefficient 



t 

m po 




(3-4) 



9 



b. The Melting Time. 



Forming the integration Eq.(3-2), by neglecting the radiation term and introducing the 
boundary condition of ^ = 1 when t = t^, yields ; 

-TJ (t - 1.) = - ! 3),p„. - p„,) (5’ - 1) (3-5) 

End of melting occurs when 0 = 0 and at that time t = t^ 



^ = 



+ 



J ( Pm -PmJ / ( 3X ( T - T„ )] 
Msf po ' ^Ms ' c' r Mm'*' 



(3-6) 



Using equation (3-4) for t^, we will find a general expression for the minimum complete 
melting time t^ 



XJr 

f po 




Msf 



MT.-tj] 



(3-7) 



The initial temperature T. of the particle, at its initial trajectory point, is determined by 
conservation of energy. We will assume that T. = 500 K, and tj^/r^^ will be found by; 

0^= (932. -500.) / (11766.- 500.) = 0.0383 

and inserting aluminum metal properties into Eq. (3-7), yields: 



t, /r = 0.0874 - 2O.23*fn(1-0 ) 

f PO ' m' 



(3-8) 
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By introducing 0^= 0.0383, into Eq. (3-8) we determine equal 66.2 sec/m . Let us 
assume a traveling time of 200 microseconds for a particle to reach the target. Now we can 
find, by substituting t^ =200 microseconds, that just 6 micron particles or smaller will 
completely melt and particles with larger diameter will have no chance to reach the igni- 
tion point. 
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4. EXPERIMENTAL RESULTS IN WATER VAPOR ENVIRONMENT 



As mentioned in chapter 2, combustion behaviour varies with the change of the envi- 
ronmental composition. Generally, the burning rate depends on the oxidizer concentration 
( Refs. 7,16 ), therefore combustion in water vapor is faster than in air at the same static 
environmental pressure and temperature. On the other hand, ignition in water vapor is 
more difficult than in air. The oxide layer deposit in wet environments is more protective 
than in dry environments and therefore ignition temperature in water vapor is 1700 K ( Ref. 
8), whereas in dry air, about 1000 K. The predominant reaction of aluminum in water has 
the form of: 



2 Al + 3 HzO = AI 2 O 3 + 3 H 2 



Wilson et al. ( Ref. 11 ) found the reaction rate constants for the above reaction. 



K = 4.916 10* Exp[ -73.5/(1.99T)] (4-1) 



Ref. 11 also mentioned that the condensed oxide smoke is carried away faster in moist 
steam than in dry air. Ref. 7 found that hydrogen attacks the oxide deposition layer causing 
porosity and, therefore, greater vapor metal fuel consumption. 

Experiments with exploding wire in steam were carried out by Kol et al., Refs. 2-4. Fig. 
3, from Ref. 4, shows an example of magnesium exploding wire and its different main com- 
bustion regions. The white brightest sphere at the center is the plasma core. Its tempera- 
tures high above 6000 K, and it contains partially ionized metal vapor. The outer spherical 
shell consists of metal vapor, where a temperature higher than the OBT is maintained. On 
the outer surface, the metal vapors meet the surrounding steam, and vapor phase reac- 
tions occur. 



Since radiative and conductive heat fluxes from the plasma core to the metal vapor shell 
are high, metal vapor temperature is expected to be much higher than the OBT. Because of 




Fig. 3 Magnesium exploding wire - region of combustion in steam 
(from Ref. 4) 

these enormous heat fluxes and high temperatures, the system can not be considered as 
an adiabatic system, and, hence, local oxide product may still appear on the gas phase. 
Results of measurements presented in Ref. 4 show that up to a Normalized Distance (ND) 
of 120 from the plasma core center, the cloud and its streaks are blue and the temperature 
is above the OBT. Beyond this distance, color is changing and temperatures, detected by 
a two color photo-pyrometry technique, indicate condensation of the oxide gas species. 
Overall maximum resolution of the measurement device is 2 micron, and, therefore, it 
could not resolve the small condensed oxide particles in the streak. As a consequence of 
the above measurement results, the author concludes that the cloud around the exploding 
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wire consists of plasma, gas phase reactant and gas products. Due to high temperatures 

in that region, metal vapor and oxide product condensation is not possible in the corona 

sphere close to the exploding center. Oxide product condensation may take place just 

after expanding away from the explosion center. As a result of the wire explosion, some of 

its molten parts are blowing-off as metal droplets and appear beyond the corona's streaks. 

Fig. 4, which presents a schematic model for emission from the metal droplet, is based on 

a color slide picture taken in a magnesium wire explosion test. In some cases, one can 

assume that hot droplets emitted from the hottest melted wire may ignite, but most of them 

do not. They simply cool and condense along their trajectory in steam. Ref. 2 describes 

2 

measurement results of the high particle deceleration (of the order of 100000. m/s ), 
caused by growth of liquid alumina on the surface of moving burning metal droplets. In 
most deceleration measurements, in Ref. 2, the solid oxide particles were the product of 
abrupt solidification of molten super cooled droplets. 

The various cases, described above, present the complexity of the phenomena of 
metal combustion under various environmental and initial conditions. In order to have bet- 
ter estimates, there is a need for a calculating tool. The next chapter is devoted to 
describing the transient combustion. 
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Fig. 4 Regions surrounding a metal exploding wire 



5. DISCUSSION OF PARAMETERS FOR MODELING THE METAL COMBUSTION 

Exploding wire and underwater shaped charge metal jet are just two instances of a 
very fast, non-steady phenomena. Even metal particles ignition and combustion in a rocket 
motor environment cannot always be considered as a quasi-steady process. A precise 
prediction of ignition time, temperatures and species existing in the various combustion 
regions can be done only if transient effects are included. 
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All quasi-steady combustion models (Refs. 6,8,12 ) assume Lewis Number unity, which 
physically means that the evaporating metal mass flux is independent of the gas phase 
reaction rate. That assumption may not be valid in any of the following circumstances: 
Reactions are kinetically controlled , strong mass convection effect and dissociation caused 
by high temperature. Those chemical and flow conditions may certainly be present at the 
beginning and the end of the combustion process of a moving metal droplet. 

The fact that the moving burning metal droplet accumulates part of its oxide product 
induces a dramatic changes in its velocity, Ref. 2. A physical model will be required to 
describe the evolution in time of the various parameters. In order to achieve that model, 
solution of the flow field, mass and heat transfer equations is needed. It is also important 
to estimate the radiation and the effect of deposition of smoke oxide cloud particles on the 
burning metal droplet. It is necessary to account for the flow as two-phase flow together 
with radiation terms. Radiation fluxes from the condensed oxide causes heat losses from 
the reaction zone and heats the burning surface. As a consequence of two-phase flow heat 
and momentum are delivered from one phase to the other, and convection terms between 
two phases change temperatures and velocities. Therefore, these effects should not be 
ignored. 

For the same kind of metal and environmental conditions, differences in results 
between the various 1-D, quasi-steady models ( Refs. 6,8,12 ) and unsteady, 2-D two-phase 
flow model are expected. The validity of including the volatile and non-volatile metals in 
one model has not yet been demonstrated. The idea of joining them into one model is 
based under the following thermodynamic facts. The difference between these two types 
of metal is mainly in the reaction zone location and behavior of oxide product. By solving 
the coupled conservation equations, one can find the distribution of enthalpy and vapor 
concentration distribution in time and space. But metal surface temperature and location 
are yet not known. By writing the mixture total enthalpy expression and including terms of 
species concentration, heat of formation, metal phase-change enthalpy, and by employing 
the parameter of vapor metal quality and the thermodynamic relationship of wall tempera- 
ture and metal-vapor concentration, one can eliminate the instantaneous surface location 
and its temperature. In the case of non-volatile metal combustion, the thermodynamic 
dependence of surface metal temperature on vapor concentration along with solution of 
conservation equations will result in the highest temperature being on the burning metal 
surface. Whereas in the case of volatile metal burning, the point of highest temperature 
may be anywhere in the field. Oxide condensation takes place whenever the local gas 
temperature is lower than OBT. Since MBT of a non-volatile metal is greater than OBT, 
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oxide condensation may appear at a large stand off from the metal surface. In contrast, for 
volatile metal, condensation of oxide may be on the metal surface. The above discussion 
shows that by applying the thermodynamic phenomena description and the conservation 
equation, a general transient model can be formulated for volatile as well as for non-vola- 
tile metal combustion. The next step will describe the governing equations of such a 
model. 
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6. GOVERNING EQUATIONS 



The various cases of exploding wire or droplet metal combustion in steam can be 
solved by one general model, which includes the transient, two-dimensional, asymmetrical 
effects. For the sake of simplicity, since we are interested in investigating the transient 
two phase flow and radiation effects, we will first consider the problem in transient one-di- 
mensional, spherical configuration. 

6.1 Energy Equation for Gas and Liquid Metal 

The equation will be written for enthalpy and will include in one region (subscript 1) 
the various gas species and liquid metal as well. Properties at each point will be pre- 
sented as dependent on species mass fraction and local temperature. Assuming constant 
static pressure in the entire region, the equation in the radial direction will be given as fol- 
lows: 



a(p,h,)/a = r*''^[r'K,(5T,/5r)]/a + o", - p,u,^(h,/^r) 



( 6 - 1 ) 



Where the mixture enthalpy is given by: 



hi = Eh.c. 
1 1 



( 6 - 2 ) 



Cj is the mass fraction of species i, and hj is given by: 



@ ' =Metal,H20,H2 



(6-3) 
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Metal enthalpy, , will account for quality of metal vapor, x, and phase change 
enthalpy,h„^g 






(6-4) 



ff n 

The heat source, Q i is the sum of other various sources, Q i. : 

» ^ 

a. Heat released by metal combustion : 



Q„ = w„H^^ (6-5) 



b. Heat released by condensation and/or solidification of oxide product : 



Q 12 



w h 



OX 



oxfg 



@T = 



T 



oxb 



(6-6a) 



Q 13 = w h 

ox 



oxfs 



@T 



T 

oxm 



(6-6b) 



Two assumptions are made here. The first is based on the idea that oxide phase 
change will occur if local temperature is lower then the OBT or lower then Oxide Melting 
Temperature (OMT). The other assumption is based on the fact that mass transfer between 
phases is a comparatively fast process. 
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c. Heat transfer from one phase to the other. 



Q 15 = ( h^A^/V )(T,-T2) (6-7) 



d. Radiation heat source due to cloud particles, 
q ^ and q are the heat radiation sources, located at the cloud and at the wall respec- 

C 

tively. Later-on, we will develop those two radiation terms. In general, q ^ stands for radi- 
ation flux, and and the radiation source will be: 



Q 14 = q /I (6-8) 



The solution of the energy equation together with the conservation equations will 
produce an enthalpy distribution for each time step. Temperature and metal metal vapor 
distributions can be eliminated from these transport equation results. Going from the far 
outer gas edge toward the metal surface, whose location has yet not been determined, the 
following sequencing logic should be adopted. For each nodal point, we need to determine 
whether liquid metal existed there at the previous time step. It is assumed that once the 
gas phase appears it will remain gas and liquid metal cannot appear again at that nodal 
point. If gas phase existed, x will be remain equal to unity, and the temperature will be 
computed from the calculated local enthalpy and gas species mass fraction. This assump- 
tion is justified by the fact that in most combustion stages, metal vapor is totally consumed 
in the reaction zone, and its concentration during the ignition period is very small. 

By assuming a constant value of C^, enthalpy for each species can be defined as fol- 
lows : 

For fully melted metal ingredient at Tj^ temperature. 
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(6-9) 



^m1 



c ( T„ - 298. ) + - 

' Mm ’ Msf 



^Ms 



For fully vaporized metal ingredient at 



*^m 3 "'"wb ' "'"nm ^ ^ ^Mfg 



( 6 - 10 ) 



Enthalpy of vapor metal ingredient when T is greater than Tj^ is given by; 



^Mg ^m3 




( 6 - 11 ) 



Enthalpy of steam ingredient at the saturation temperature will Include its heat of for- 
0 

mation, h^ , at 298. K and phase change latent heat, h^^^. 



h = 

wsv 



+ C„ (T^b-298-) + h 






wfg 



( 6 - 12 ) 



Water vapor ingredient enthalpy at temperature above the saturation point will include 
h . 

wsv 
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(6-13) 



h = 



WV 




Hydrogen product of combustion will always appear in the gas phase, and its enthalpy 
term will be ; 



= "p ( T - 298. ) 

n 



(6-14) 



According to Eq. 6-2 the total enthalpy of the gaseous mixture is defined as ; 



h = h„c„ 



h c -t- 

WV WV 






(6-15) 



After introducing the terms for each ingredient enthalpy Eqs. 6-11,13,14 and defining 




CmC,j+C c -hcc„ 

pMg M p^ WV Pj^ H 



(6-16) 



The temperature at the nodal gas phase can be eliminated, and its expression is : 



T = [''■( - c )c„ - (h_ - 373.)c^ • c 



'^Mg 



c„298.] / ( Cp ) (6-17) 

g 
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If at the previous time step, x was less than 1, or even equal to 0., the procedure for 
determining a new value for x, should be as follows : 

By knowing the local metal vapor mass fraction and the total local molecular weight, 
assume the local temperature to be equal to 



For aluminum A= 36560.8 and B= 25.016. 

Metal vapor quality value, x, can be eliminated from the total enthalpy expression. By 
assuming that the nodal ingredient mixture consists of saturated metal vapor, water vapor 
and hydrogen, and assign the local temperature equal to x remains the only unknown 
and, therefore, can be eliminated: 



If the resulting value of x, based on those assumptions, is larger than one, then the local 
temperature should be greater than the assumed temperature and a new tentative 
value of temperature should be determined by using Eq. 6-17. If , at the first iteration x 
becomes equal to or greater then zero, then the temperature is set equal to Tj^ and the 
surface location is determined. Further inside the metal, temperatures are determining by 
setting x equal to zero. 



= A/[B-£n((c„M^)/M„)] 



(6-18) 




(6-19) 



T = (h - h 1 + c T ) / c 



( 6 - 20 ) 
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The thermodynamic relationship for the metal vapor density shows that increasing of x will 
decrease the density Eq. 6-21. Metal surface vaporization cause a local density 
decrease and since the continuity law must also be conserved, the net metal vapor will 
expand outward to the gas surrounding. 



Pm = pMg / n + (1 - X )( - 1 )] 



( 6 - 21 ) 



6.2 Determining the Effective Radiation Heat Source 

Considerabie amounts of radiating heat fluxes govern the metal combustion process in 
the case of large volatile particles. In most other cases the surrounding gas can be con- 
sidered as an opticaliy thin medium. When oxide particies appear, the media is not trans- 
parent any more. The number of particles per unit voiume N, the particie size r^, specific 
particle emissivity, cloud thickness I, and the relative geometric location of the particle 
cloud determines the radiation fluxes. 

Intensity attenuation of radiation propagating through a particle cloud is given in Ref. 
13 as: 



i/io = EXP( -Gi ) 



( 6 - 22 ) 



Where the extinction coefficient g is a function of N and r . Measurement results described 

p 

in Ref. 13 demonstrate the validitiy of the relationship between these parameters. We 
present in Fig. 5 some data results of aluminum oxide cloud, which will be used as a 
data base for our model. 
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Particle radius 




Fig. 5 Correlation of all spectal data. ( from Ref, 13 ) 

Solution of the gas phase chemical reaction along with the other conservation equa- 
tions will produce the oxide product distribution in the field. By assuming an average par- 
ticle size rp, and knowing the local cell density, we are able to determine the number of 
local oxide particles in that particular volume cell ; 



N. = (3./4.;trp') (p./pj = ^2V/(4rtrp&S"\/3) 



(6-23) 
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By accounting for all nodal N.'s and their radial position, it is simple to find the total parti- 
cles cloud depth and its average number of particle per unit volume. 



N = ( LN.AV. ) / ( EAV.) 

\ 11'' l' 



(6-24) 



Using the chart in Fig. 5, we can find the value of e then the overall cloud emissivity is 
determined from : 

Ep = 1 - EXP ( - gI ) (6-25) 

Fig. 6 shows directions and radiation fluxes interchanges between the cloud, the metal 
surface and the environmental gas. For practical purposes it is assumed that all media 
involved are grey. As the gas is optically thin, and particles are mainly concentrated in a 
narrow zone (usually the flame zone), isothermal radiating surfaces is a resonable 
assumption. Total particle cloud radiation is diffused toward the metal surface and to the 
outer gas space. Part of the radiation emitted from the cloud and impiging on the metal 
surface, will be reflected back to the cloud by the grey metal wall and absorbed either by 
the particles in the cioud or by the gas. The net result of radiation fluxes through the par- 
ticles cloud is also a function of the mean gas and particle cloud temperatures, T^, T^^ and 
wall metal temperature T^ 



(q /<r) = aET F s.+oeT s 

' p w w w-c 1 p g g o 



-1- a (e T ^(1-a )F s. -e T (6-26) 

p' P OX ' w' W“C 1 p ox ' ' 



The s. and s^ are parameters of surface ratios, as shown in Fig. 7 and for spherical geom- 
etry they are given by : 



s. 

1 



= A . / ( A . -F 



A ) = 

CO ' 



/ (r . + r ) 

' o 1 no • 



(6-27) 
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( 6 - 28 ) 



A / ( A . 

CO ' Cl 



A ) = 

CO ' 



/(r . + r ) 



The geometrical view factor is given as: 



F 

c-w 



[1./ (E 



.1 



(V' - 1 V Aei))] 



( 6 - 29 ) 



The metal surface absorbs radiation from the particle cloud and the gas while emitting its 
radiation heat flux outward. 



(q. 



net 



IfS)- 



a € T 

W P ox 



c-w 



a £ T 

w g g 



-E T 

w w 



( 6 - 30 ) 



Those two radiating fluxes divided by the appropriate characteristic length, which is the 
particle cloud thickness, provide the local corresponding heat source that should be taken 
into account as source terms in the energy equation. 



27 



gas 




^ Ul I u 



Fig. 6 Radiation heat fluxes exchange between the oxide cloud 
and the nnetal surface. 

6.3 Energy Equation for the Oxide Particle Phase 

The smoke oxide product of a volatile metal combustion requires an energy equation, 
as it is a second phase in the field. Since the small oxide particles are dispersed in the 
gas, conduction heat exchange between them can be neglected. The heat transfer time 
constant for the largest oxide particle (which is less than 2 micron), in a laminar flow is 
less than 10 microseconds, therefore the radiation heat losses of the particle cloud may be 
incorporated as a source term in the gas phase energy equation. If local temperature is 
smaller than the OBT, then' immediate condensation of all the new amount of the oxide 
product will occur. 
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Fig. 7 Geometrical relationship between the cloud and 
metal droplet. 

If ^2 represents the volume fraction of phase 2, then the contact surface 
the oxide particles with an average radius r^ in a volume V is given by Ref. 14 as: 



A = 3^2V / r 

S P 



between 



(6-31) 
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The source term for that phase will include internal-phase heat transfer and release of 
latent heat. Ref. 14 presents the source term as : 



Qz" = X (A /V)( T 2 - T, ) - w "h , 

c' s ' OX oxfg 



(6-32) 



The energy equation will, therefore, be presented as ; 



^(P 2 h 2 )/^ = Q 2 - p2U2^h2/5r 



(6-33) 



The oxide temperature can be directly found from its enthalpy h 2 



T = (h 2 -h -c 298.) / c 

ox ' ox P ' F 



(6-33a) 



6.4 Mass Transfer Equation for Each Species 

Analytical descriptions of species mass fraction and enthalpy transport are of the same 
kind as the general diffusion equation. At some time, when Lewis Number equals 1., their 
non-dimensional solutions are the same. The distribution in time and space of both vari- 
ables, Cj and h, are given by the solution of a diffusion equation. With the assumptions of 
quasi-steady combustion and unit Lewis Number unity, the numerical problem reduces to 
an algebraic set of 9 equations, (Ref. 8). But in many particular combustion cases, such as 
an exploding wire or a composite propellant, high temperatures cause dissociation, so that 
assumption is not valid any more. Simultaneous solution of the diffusion equation is 
needed with careful attention to various source terms and temperature dependence of the 
species properties. The source term in the mass transfer equation is the mass production 
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or consumption of the each particular species. At each time step and location, the source 
term is determined by local thermodynamic condition. If the local temperature is not high, 
enough the species source term will be determined by the reaction rate as: 



u 0 

= B T"EXP(-E^RT)^m wv 



(6-34) 



If the reaction rate consumes more than the species mass flux into the local volume cell, 
then the local consumption and production are controlled by the 'fuel' and 'oxidizer' diffu- 
sion rates. If the local mixture ratio f. Eg. 6-35a is less than the stochiometric mixture ratio 
fgj (Eq. 35b), then the new value of the local oxidizer, (water vapor), mass fraction, will be 
governed by the mixture fraction and all the fuel, (metal vapor), will be consumed. A new 
value for fuel mass fraction will be set and complete consumption of the oxidizer will occur 
when f is greater than f^^. This procedure is presented graphically in Fig. 8, and Eqs. 36 
show its analytical representation. 

f = 1./(1+c^/c„) (6-35a) 

f^^ = 1./(1. + (cjc^)j (6-35b) 



For 0<f%ofg^ 




For 1 ftf^, 



c = 1.- f/f . 

wv st 

c =0. 

wv 

c„=(f-f,,)/(1-f,,) 



(6-36a) 



(6-36b) 



The products mass fractions are determined according to the stochiometric ratios i. One 
Kgm of aluminum will produce i^^= 1.88653 kgm of alumina and 1^^ = 0.11898 kgm of hydro- 
gen, in the process it will consume i^= 0.998 kgm of water. The net rate of oxide produc- 
tion will depend on the mixture condition and mass fraction of the reactants. 

For a fuel rich mixture the oxide mass fraction source is: 
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Fig. 8 Mass fraction distribution with mixture fraction. 

= (‘ox / 



(6-37a) 



and, for the fuel lean case, it is: 



// 

w 



ox 






(6-37b) 
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Since the hydrogen product will always appear in the gas phase, we are allowed to set its 
mass fraction equal to: 




(6-38) 



The diffusion equation for the mass transport of species j is given by; 



5(pic^)/^ = r'^^[r^pj^Dj^(^c./5r)]/5r -h w^ - pjUj(^c75r) 



(6-39) 



The density is the density of the local mixture, and 0.^ is the species diffusion to the 
mixture. Late in chapter 6.7, we will show how to specify the local properties. 

6.5 The momentum equations for each phase 

The two phases have the same type of the momentum equation with a source term 
that stands for the momentum exchange between them. We assume the second phase 
consists of solely oxide particles. Since the particles are dispersed in the gas, they do not 
exert friction forces on each other, and, therefore, we may ignore the momentum transport 
within that phase. We now define the source term of momentum interphase transfer. From 
Stock's flow around one particle we know that the friction force, ,is given by : 



F = 6niir ( Ui - U 2 ) 

P P 



(6-40) 



Now we define F as the force F per unit 



p p 
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particle volume. 



F '' = F /(4rtf "/3) = 4.5 \i!r " 

p p ' p ' p 



(6-41) 



If ^1 and ^2 are the volume fractions of phase 1 and 2 respectively, then the momentum 
interphase exchange per unit volume and time will be presented as follows: 



IVIi2 — - IVI 21 “ F V(ui*U2) 

P 



(6-42) 



The radial one-dimensional momentum equation for each phase will then be given by: 



pi(^Ui/^) = r''^^ [r^p,p,(^u,/^r)]/^r - p,u,(^u,/^r) -I- M 21 (6-43) 

P2(^U2/^) = - P2U2(^U2/^r) + Mi2 (6-44) 



Since we solve the gas and the metal as one phase by ignoring surface tension forces 
between them, the metai droplet may be smeared into the gas. In order to prevent that, 
both phase velocities inside the metal should be set equal to zero. For small enough time 
step the procedure of freezing velocities will capture the oxide on the metal surface, and 
therefore simulate deposition. 
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The total diffusion of species j into a mixture which contains a mole fraction x of that spec- 
ies, is defined by 



- ( 1 - Zj ) / (!:(Zk ! 



(6-49) 



The mole fraction x-. can be found from c. , M. and M. as follows: 



X- = c.M. / M. 

^3 3 -t 3 



(6-50) 



The mixture conductivity and viscosity coefficients can be found from gas tables Ref.15. 
They are expressed as a fourth degree temperature polynomial for each species present in 
the mixture. The conductivity k. , the viscosity and heat capacity c will be presented 

3 3 Pj 

as: 



k.= Ek^^T^ 



(6-51 a) 
(6-51 b) 




Ec T^ 
Pl3 



(6-51 c) 
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In order to evaluate the mixture viscosity, a bi-molecular viscosity parameter SI., should be 
defined : 






( 6 - 52 ) 



H is defined as: 



H = S ^ 



( 6 - 53 ) 



The mixture conductivity will be written as: 



K=Zfaikj/(i:(x,n.,))i 



( 6 - 54 ) 



The density pi is defined by the ideal gas equation of state. 



P. = p/(RTi) 



( 6 - 55 ) 



c was defined previously, in 6-1 when discussed the gas energy equation. The next step 

^Mg 

should be to cast the model, presented above, in the form of a computer program. 
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7. SUMMARY 



Preliminary calculation for aluminum particles from a shaped charge revealed that the 
maximum size of a particle that may reach ignition conditions while moving at speed of 
7000. m/s during 200 microseconds in a steam environment, is just 6 micron. 

Exploding wire experiments in water vapor chamber is an adequate technique for find- 
ing the most combustion effective metal liner of underwater shaped charge. However, it 
does not simulate the phenomena of ignition and combustion of high speed moving metal 
particles in the U/W S/C. The need for precise prediction in any of these cases requires a 
solution, which will account for the transient effects along with two phase flow and radia- 
tion fluxes emitted from the oxide cloud. The combustion process parameters that will 
generalize the model for both the volatile and non-volatile metals cases are: the metal 
vapor quality, surface temperature and metal vapor surface concentration. By writing the 
energy equation for total enthalpy that include those parameters and accounting for local 
metal and oxide phase change temperatures, the solution will be valid for both type of 
metals. The conservation equations were written for the entire region, gas and liquid. The 
technique of iterating in order to find the metal surface was described. Radiation is 
accounted for as source terms in the oxide and at the metal surface. 

Such a solution for a general type of metal burning will simplify the computer code 
program for the various cases. 
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